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The three and five-dimensional convex sets of two-level complex and quaternionic quantum sys- 
tems are studied in the Bayesian thermostatistical framework introduced by Lavenda. Associated 
with a given parameterization of each such set is a quantum Fisher (Helstrom) information ma- 
trix. The square root of its determinant — adopting an ansatz of Harold Jeffreys — provides a 
reparameterization-invariant prior measure over the set. Both such measures can be properly nor- 
malized and their univariate marginal probability distributions — which serve as structure functions 
— obtained. Gibbs (posterior) probability distributions can then be found, using Poisson's integral 
representation of the modified spherical Bessel functions. The square roots of the (classical) Fisher 
information of these Gibbs distributions yield (unnormalized) priors over the inverse temperature 
parameters. 

PACS Numbers 05.30.Ch, 03.65.-w, 02.50.-r 

In this letter, we apply the Bayesian thermostatistical framework of Lavenda Jl|-^| to the two- level quantum systems, 
both complex and quaternionic Q in nature. A two-level quaternionic system can be represented by a 2 x 2 density 
matrix 



1 / 1 + z x — iy — ju — kv 

2 \x + iy + ju + kv 1 — z 



(1) 



where i 2 — j = k 2 = —1, ij = —ji — k, jk = —kj = i and ki = —ik = j. Setting u = v — 0, we obtain the familiar 
(Pauli matrix) representation of the two-level complex systems sec. 4.2], in which the points (x, y, z) lie within the 
unit ball (x 2 + y 2 + z 2 < 1). 

In the complex case (u = v = 0), one can — using the concept of a symmetrized logarithmic derivative []6|-|ll"f 
associate with the given parameterization, the quantum Fisher information matrix, 



(1 — x 2 — y 2 — z 2 ) 



y 2 — z 2 xy xz 

xy 1 - x 2 — z 2 yz | . (2) 

xz yz 1 — x 2 — y 2 



In the quaternionic instance, employing the relations between the Pauli matrices and the quaternions Q p. 495] 
[ p2| , p. 197] to generate a complex 4x4 complex density matrix and then finding the corresponding symmetrized 
logarithmic derivatives |ll| — one has 



(1 — u 2 — v 2 — x 2 — y 2 — z 2 ) 



/ g{{u}) uv ux uy uz \ 

uv g{{v}) vx vy vz 

ux vx g({x}) xy xz 

uy vy xy g({y}) yz 

\ uz vz xz yz g({z}) I 



(3) 



where g({u}) denotes 1 — v 2 — x 2 — y 2 — z 2 , ... (The inverses of (||) and (||) provide — by the Cramer- Rao inequality 
- lower bounds Jin the sense of negative defmiteness] on the covariance matrices of the parameters @j@].) The 
determinants of (|^) and (||) are, respectively, 

(4) 



(1 — x 2 — y 2 — z 2 ) 
and 



(1 — u 2 — v 2 — x 2 — y 2 — z 2 ) 
(These are inversely proportional to the determinants obtained from ([!]) itself.) 



(5) 
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In the classical/nonquantum situation — in which probability distributions rather than density matrices are studied 
— Harold Jeffreys proposed that one employ the square root of the determinant of a Fisher information matrix as 
a (reparameterization-invariant) prior measure (volume element) over a family (Riemannian manifold) of probability 
distributions ]l3]-[l5| . Adopting this ansatz to the quantum mechanical case, one can — normalizing the square roots 
of (^) and (^|) over the three and five-dimensional unit balls — derive the prior probability distributions, 

' {x 2 + y 2 + z 2 <l) (6) 



and 



-y* - z 2 )V^ 
2 

7T 3 (1 — u 2 — v 2 — x 2 — y 2 — z 2 ) 1 / 2 



(u 2 + v 2 + x 2 + y 2 + z 2 < 1) 



(7) 



The three bivariate marginal probability distributions of (g) are uniform over unit disks, while the five quadrivariate 
marginal probability distributions of (0) are uniform over unit balls in four-space. Further integrating over all but 
the last coordinate, one arrives at a univariate probability distribution, having the form, in the complex case, 



and, in the quaternionic instance, 



2(1 -z 2 ) 1 / 2 



:i-z 2 ) 3 / 2 

3tt 



(-1 < z < 1) 



(-1 < z< 1) 



(8) 



(9) 



We consider (|J) and (j|) to be (normalized) structure functions in the sense of Lavenda fl|-f| . 
Since (||) and (||) are proportional to expressions of the type, 

(1 _ 2 2 r -i/ 2) (n = 1)2 ) (10) 

one can employ Poisson's integral representation of the modified spherical Bessel functions (cylinder functions of half 
integral order) |Ig|] , 



/»(/?) 



(f) r 



exp 



{-(3z)(l-z 2 ) n - 1/2 dz 



(11) 



V¥r(« + i) j_! 

to introduce thermodynamic considerations |]l7|-[l9|l. One then has (properly normalized) Gibbs distributions of the 
form, 

expRfeXfHl-z 2 )- 1 / 2 

I„(/3)^r(n+i) ' 1 ' 

where fi serves as the inverse temperature parameter. For (3 = — 1, Fig. |l| shows these distributions for n = 1 (complex) 
and n = 2 (quaternionic), the latter curve having the higher peak. 
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FIG. 1. Gibbs distributions dia) for f3 — —1. The higher-peaked curve corresponds to the quaternionic (n = 2) case. 
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Fig. g displays the analogous results for (3 = 5, with the curve for the complex case now having the higher peak. 

probability density 
2.5- 




-1 -0.5 0.5 1 

FIG. 2. Gibbs distributions @ for f3 = 5. The lower-peaked curve corresponds to the quaternionic (n = 2) case. 

In Fig. H and || are shown, respectively, the expected value of z (< z >) and the variance about < z > as a function 
of /3, with the quaternionic curves being the flatter ones in the two graphs. 
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FIG. 3. Expected value of z as a function of (3. The quaternionic (n = 2) curve is the flatter one. 
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variance 
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FIG. 4. The variance about the expected value of z as a function of (5. The quaternionic (n = 2) curve is the flatter one. 

In Fig. [5] are displayed the relative entropies of the Gibbs distributions (12) with respect to the uniform distribution 
(i) over z £ [— 1, 1], with the quaternionic curve having the greater minimum. 
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FIG. 5. Relative entropy of Gibbs distributions with respect to a uniform distribution. The quaternionic (n = 2) curve is 
the one having a greater minimum. 



Since for n = 1,2, in particular, the Gibbs distributions (|12|) form two families of probability distributions, each 
parameterized by /3, one can find the (classical) Fisher information associated with them p],p|,p^|. This is accomplished 
by computing the negatives of the expected values (relative to (|l2|)) of the second derivatives of the logarithms of ( |l2] ) 
(n = 1, 2) with respect to (3. The required integrations were performed numerically. The square roots of the obtained 
Fisher information statistics are exhibited in Fig. || for j3 £ [—10, 10], with the quaternionic (n = 2) case having the 
smaller maximum. 
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FIG. 6. Unnormalized Jeffreys' priors over inverse temperature parameter /3. The lower-peaked curve corresponds to the 
quaternionic (n = 2) case. 

These serve as the Jeffrey's (unnormalized) priors over the inverse temperature parameter (/?), in the sense of 
Lavenda's pioneering investigations (l]-|!| . (Fig. || was also reproducible, with MATHEMATIC A, using exact / numerical 
methods, but only by taking the required integrations over z G [—1,1] — which themselves proved to be problematical 
— to be the results of integrations over z € [0, 1] and the addition to them of the outcomes of substituting —/3 for (3 
in them.) 

For literature supplemental to [0j|] pertaining to the topic of temperature fluctuations, see [po|-p8| . The quantum 
Fisher information metric || employed here to generate the structure (prior) functions (^) and (|9|) has been shown 
to be simply proportional to the Bures metric — which extends to the mixed states, the Fubini-Study metric on 
the pure states of a quantum system. 
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